Density distributions of superheavy nuclei 
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We employed the Skyrme-Hartree-Fock model to investigate the density distributions and their 
dependence on nuclear shapes and isospins in the superheavy mass region. Different Skyrme forces 
were used for the calculations with a special comparison to the experimental data in 208 Pb. The 
ground-state deformations, nuclear radii, neutron skin thicknesses and a-decay energies were also 
calculated. Density distributions were discussed with the calculations of single-particle wavefunc- 
tions and shell fillings. Calculations show that deformations have considerable effects on the density 
distributions, with a detailed discussion on the 292 120 nucleus. Earlier predictions of remarkably 
low central density are not supported when deformation is allowed for. 

PACS numbers: 21.60.Jz, 21.10.Ft, 21.10.Gv, 27.90.+b 



I. INTRODUCTION 

Proton and neutron density distributions, the associ- 
ated root-mean-square (rms) radii R p and R n , and the 
neutron skin thickness AR = R n — R p provide funda- 
mental information on nuclear structure. For example, 
halo nuclei arc charactered with long tails in density dis- 
tributions. The density is a direct probe of the size of an 
atomic nucleus, and plays an important role in the cross 
sections of nuclear reactions. While charge densities can 
be measured from the elastic scattering of electrons, neu- 
tron densities are largely unknown. The parity-violating 
electron scattering has been suggested to measure neu- 
tron densities Q. With the experimental method, the 
neutron densities of more nuclei can be expected to be 
measured. Theoretically, the calculated neutron skin 
thickness can be model dependent 0j- 

Recently, Horowitz et at. |j| studied the relationship 
between the neutron skin of the spherical double magic 
nucleus 208 Pb and the properties of neutron-star crusts, 
showing the importance of the knowledge of nucleon den- 
sities in understanding the equation of state of neutron- 
rich matter and therefore the properties of neutron stars. 
The heavy nucleus 208 Pb has been measured to have a 
neutron skin thickness of about 0.15 fm Q. With increas- 
ing mass number, the neutron excess becomes larger in 
general and it's natural to think that superheavy nuclei 
provide the largest neutron excesses. The heaviest nuclei 
also have large proton numbers and thus large Coulomb 
repulsive forces that push the protons to larger radii and 
therefore change density distributions. Novel density dis- 
tributions were predicted in extraordinary A > 400 nu- 
clei that have bubbles and showed coupling effects be- 
tween density distributions and shell structures [E IE • 
The nucleon density in a bubble is reduced to be zero. 
Semi-bubble nuclei were also suggested with the consid- 



erable reduction of central densities for the Z > 120 nu- 
clei @ located around the center of the predicted island 
of stability of superheavy nuclei. Bender et al. have 
investigated the density distributions of superheavy nu- 
clei with the restriction of spherical shapes Re- 
cent progress in experiments are motivating the struc- 
ture study of superheavy nuclei [E llOj . Many theoreti- 
cal works have inve stig ated the pro per ties of superheavy 
nuclei [3 13 13 El 13 13 111 13 13 E3, such as 
shell structure, a decay and spontaneous fission. Exper- 
iments have also provided the structure information of 
superheavy nuclei by the in-beam study of spectroscopy 
[2lL |3 13 El- In the present work, we investigate 
the density distributions of superheavy nuclei and related 
structure properties, with deformation effects taken into 
account. 



II. CALCULATIONS 

The deformed Skyrme-Hartree-Fock model (SHF) [3] 
was used in the present investigation. Pairing correla- 
tions are treated in the BCS scheme using a 5-pairing 
force, V pail = -V q S{n - r 2 ) |3 0. The pairing 
strength V q (q=p, n for the protons and neutrons, re- 
spectively )has been parameterized throughout the chart 
of nuclei [23], but the actual values are dependent on 
Skyrme forces chosen. The detailed values of the pairing 
strengthes can be found in Ref. Q . 

Calculations are performed in coordinate space with 
axially symmetric shape. The density distribution of pro- 
tons or neutrons is given in the two-dimensional form as 
follows. 

p(z, r) = J2 2 «*W*(*. r )f + 1^ (*> ^l 2 ) (1) 



where, tp^ and ip k 



are the components of the wavefunc- 
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tions with intrinsic spin s z =+-| and — |, respectively, and 
v\ is the pairing occupation probability of the k-th orbit. 
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FIG. 1: The calculated proton (p p ) and neutron (p n ) density 
distributions of the spherical nucleus 208 Pb, compared with 
the experiment for the charge density |3^| . 



FIG. 2: The calculated proton and neutron densities for 
spherical nucleus 298 114. 



A. Density distributions of spherical nuclei 



The ground states of most superheavy nuclei are expected 
to have axially symmetric or spherical shape |l3j . In this 
paper, we consider the most important axially symmetric 
deformations, /?2 and (3^. 

In the present work, we investigated the densities and 
related structure problems of even-even superheavy nu- 
clei with Z=104— 120. In the SHF calculations, results 
are in general parameter dependent. For example, the 
SkI3 force predicts 292 120 for the next magic nucleus 
beyond 208 Pb, while SLy7 and SkI4 predict 310 126 and 
298 114 for the magic nucleus, respectively 0, [n|. To 
make comparison, we used the different sets of param- 
eters SLy4, SLy7 H|, SkI3 and SkI4 [29|. These sets 
of parameters have been developed recently with good 
isospin properties, and we note that they have been rec- 
ommended by Rikovska Stone et al. for their ability 
to describe realistic neutron stars and the properties of 
asymmetric nuclear matter |30j . In the superheavy mass 
region, Skyrme parameter sets can reproduce experimen- 
tal binding energies within a few MeV |31j and a-decay 
energies within a few hundred keV [TJ. Table I lists 
the properties of the ground states calculated with SkI4 
for experimentally known even-even superheavy nuclei 
and the predicted magic nucleus 298 114. The calculated 
a-decay energies agree with experimental data within a 
few hundred keV (The largest difference with data is 640 
keV in 266 Hs). Experimental binding energies can be re- 
produced within rs 4 MeV for the nuclei listed in Table 
I with the SkI4 force. The obtained neutron-skin thick- 
nesses are smaller than the calculations by the relativistic 
mean- field (RMF) 0. (It was pointed out that RMF 
calculations overestimate neutron-skin thicknesses |34j-1 



In order to test potential parameters in the calcula- 
tions of densities, we calculated the density distributions 
of the spherical doubly magic nucleus 208 Pb using the dif- 
ferent sets of parameters. For 208 Pb, the charge distribu- 
tion has been measured by electron scattering j3f|. Fig.l 
shows the calculated density distributions with compari- 
son with the experimental charge density. It can be seen 
that the proton density given by the SkI4 force is clos- 



TABLE I: The SkI4 calculations for the experimentally syn- 
thesized even-even superheavy nuclei. The calculated a- 
decay energies (Q^ 1 ) are compared with experimental values 
(Qa Xpt ). The experimental data are taken from [3rjll . The 
calculation for the predicted magic nucleus 298 114i84 is also 
listed. 



nuclei 
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0* 


Rn 


Rp 


AR 




^cxpt . 








(fm) 


(fm) 


(fm) 


MeV 


MeV 


254 Rf 


0.29 


0.06 


6.05 


5.96 


0.09 


9.35 


9.38 


256 Rf 


0.29 


0.04 


6.07 


5.97 


0.10 


9.47 


8.96 


258 Rf 


0.28 


0.03 


6.09 


5.98 


0.11 


9.07 


9.25 


258 Sg 


0.28 


0.03 


6.08 


5.99 


0.09 


10.04 


9.70 


260g g 


0.28 


0.01 


6.10 


6.00 


0.10 


9.65 


9.93 


266 Sg 


0.26 


-0.03 


6.15 


6.03 


0.12 


8.28 


8.88 


264 Hs 


0.27 


-0.02 


6.12 


6.03 


0.09 


10.24 


10.59 


266 Hs 


0.26 


-0.03 


6.14 


6.04 


0.10 


9.70 


10.34 


270 Hs 


0.25 


-0.06 


6.18 


6.06 


0.12 


8.94 


9.30 


270 110 


0.25 


-0.05 


6.16 


6.07 


0.09 


11.32 


11.20 


284 112 


0.15 


-0.07 


6.25 


6.12 


0.13 


9.65 


9.30 


288 114 


0.13 


-0.08 


6.27 


6.14 


0.13 


10.36 


9.97 


292 116 


0.06 


-0.02 


6.27 


6.15 


0.12 


10.64 


10.71 


298 114 


0.0 


0.0 


6.33 


6.16 


0.16 


7.4 
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FIG. 3: The SkI4 calculated \ipk\ 2 distributions of the proton 
lii3/2, lhg/2 and 2f 7 / 2 orbits in the case of spherical 298 114 
nucleus. The dashed line stands for the summed contribution 
of the orbits. 

est to the experimental measurement. The rms radii of 
208 Pb are calculated with the SkI4 force to be 5.43 fm 
for the protons and 5.61 fm for the neutrons, leading 
to a neutron-skin thickness of 0.18 fm which is slightly 
larger than the values of 0.16 fm given by the SLy4 and 
SLy7 forces. These results agree with the 0.15±0.02 fm 
from the recent antiprotonic atom experiment 0. The 
SkI3 force gives a larger neutron skin thickness of 0.23 fm 
compared to the other three forces, which may be due to 
the similar behavior of SkI3 force to the RMF model Q . 
It needs to be pointed out that nuclear ground-state cor- 
relations (see, e.g., |33, H3) can have visible effects on 
nuclear properties, such as energies and densities. The 
oscillations observed in the calculated density distribu- 
tions in Fig.l could be reduced when the correlation is 
taken into account. Such correlations, which go beyond 
the mean-field approximation, are not included in the 
present work. 

Fig. 2 shows the calculated densities of the spherical 
nucleus 298 114. This nucleus was predicted to be the 
next doubly closed shell nucleus by a Macro-microscopic 
model [H and SHF with SkI4 force 0. For 298 114, the 
calculations with SkI4, SLy7 and SLy4 give similar den- 
sity distributions. It can be seen that the charge density 
distribution of 298 114 has a central depression. The cen- 
tral density depression has been predicted to exist widely 
in the spherical superheavy nuclei ||. Fig. 3 displays 
the calculated square wavefunctions of the proton H13/2, 
lhg/2 and 2f7/ 2 orbits that locate the 82<Z^114 closed 
shell. It can be seen that the high-j orbits have den- 
sity contributions in the nuclear surface region. This is 
consistent with the classical picture in which orbits with 
large angular momentum locate at surface. In Fig. 4, we 
show the proton densities of 208 Pb and 298 114 for com- 
parison. The proton density of 298 114 is decomposed into 
two parts: i)the contribution from proton orbits below 
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FIG. 4: The SkI4 calculated proton densities of 298 114 and 
208 Pb. In the 298 114 nucleus, the dashed line presents the con- 
tribution from the orbits below the Z=82 closed shell, while 
the dot-dashed line presents the contribution from the high-j 
orbits of lii3/2, lhg/2 and 2f 7 /2 that construct the next closed 
shell. 



the Z=82 closed shell; and ii) from the orbits in the next 
closed shell with 82<Z<114. The contribution from the 
orbits below Z=82 has a similar behavior to the charge 
distribution of 208 Pb, without central depression. The 
82<Z<114 orbits (that have high-j values) have contri- 
butions in the surface region of the nucleus, leading to a 
central depression in the proton density of 298 114. 

For the neutrons at spherical case, the high-j orbits 
of 2g 9 / 2 , lin/2, lji5/2 and 2g 7/2 occur in the region of 
N=126— 172. The low-j orbits of 4s!/ 2 and 3d 3 / 2 are at 
N=178— 184. Fig. 5 shows the density distributions for 
N=172— 196 and Z=114 with assumed spherical shape. 
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FIG. 5: The spherical density profiles of the nuclei with 
Z=114 and N=172-196 calculated with the SkI4 force. 
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Indeed, these nuclei were predicted to be nearly spherical 
in their ground states [IT| It can be seen that proton 
densities have central depressions and neutron densities 
become centrally depressed for N<178. 



B. Densities of deformed superheavy nuclei 

The shell fillings of nucleons are sensitive to the defor- 
mations of nuclei. Therefore, density distributions should 
be expected to be shape dependent. Most superheavy nu- 
clei known experimentally are believed to have deformed 
shapes. Fig. 6 displays the calculated density distribu- 
tions for the Z=110 isotopes with N=160, 170 and 180 
nuclei. The equilibrium deformations are determined by 
minimizing calculated energies. The deformations deter- 
mined with the SkI4 force are /3 2 =0.25, 0.19 and 0.05 
£ or 270,280,290 2 respectively. To check the possible pa- 
rameter dependence, we also used the SLv7 force 1 to cal- 
culate the densities, shown in Fig. 6. The deformations 
determined with SLy7 are /3 2 =0.25, 0.17 and 0.05 for 
270,280,290-qq^ reS pectively. In order to sec deformation 
effects, we calculated the densities assuming spherical 
shape. It can be seen for 270 > 280 no that the densities 
become more centrally depressed in spherical cases. We 
also see some difference between the densities along the z- 
and r-axes. Fig. 7 shows the two-dimensional proton den- 
sity distribution for the 280 110 nucleus. The density has 
two humps in the z-axis. The double-hump distribution 
has been suggested experimentally, e.g., in the deformed 
166 Er and 176 Yb [37|- Our calculations with SkI4 and 
SLy7 show that such doublc-hump phenomenon is rela- 
tively pronounced for nuclei around 280 110. The macro- 
microscopic calculations show that the even-even nuclei 
around N=170 are particularly unstable against sponta- 
neous fission (see Fig. 9 in This could be related to 
the double-hump distributions in these deformed nuclei. 

Fig. 8 shows the distributions of the square wavefunc- 
tions of the j z = 1/2 orbits in 280 110, modified with 
pairing occupation probabilities (see Eq.(l)). These or- 
bits are above the Z=82 and N=126 shells for the protons 
and neutrons, respectively. Shown as in Fig. 8, the high-j 
low- j z orbits have important contributions to densities in 
the surface region of the z-axis. High-j low-j z orbits have 
strong prolate-driving effect. Hence, sufficient number of 
high-j low- j z orbits occupied can result in prolate shapes 
and double-hump densities. In 166 Er and 176 Yb that were 
suggested experimentally to have double-hump densities, 
the low-j z orbits of the proton lhix/2, lg7/2 subshells are 
occupied. 

Density distributions given by different Skyrme forces 
can differ as shown in Fig. 6. The SLy7 predicts larger 
central depression in neutron densities and less central 
depression in proton densities than the SkI4 force. This 
difference also occurs in the calculation of 208 Pb, see 
Fig.l. For different parameters, the opposite behavior 
of proton and neutron distributions reduce the difference 
in the total (proton+neutron) nuclear density distribu- 



tions. The origin of the opposite behavior would be due 
to the self-consistent coupling between protons and neu- 
trons in the SHF model, to approach the nuclear density 
saturation j^. The SkI4 force is better in reproducing 
the density of 208 Pb, compared to other Skyrme forces 
used in the present investigation. However, the SkI4 force 
was pointed out to have larger spin-orbit splittings in the 
calculations of single-particle level schemes for the super- 
heavy region Q . 

C. The 292 1 20 nucleus 

292 120 is an interesting nucleus that was predicted to 
be a doubly magic Q and spherical semi-bubble nucleus 
0. We calculated the energy curve with the SLy4, SLy7, 
SkI3 and SkI4 parameters, see Fig. 9. The results of SLy4, 
SLy7 and SkI4 are close to the Hartree-Fock calculations 
with Gogny force [(|. The SkI3 force gives a shallow 
spherical minimum. The SkI4 calculation predicts three 
minima at fa =0.11, —0.12 and 0.52 (Supcrdcformations 
in superheavy nuclei have been discussed by Ren et al. 
CU D3 ) • It needs to be mentioned that the present calcu- 
lations are restricted to axially symmetric shapes without 
considering the possibility of triaxiality. The inclusion of 
the triaxial degree of freedom could alter the shallow min- 
ima. Fig. 10 shows the SkI4 calculated density distribu- 
tions at the prolate and oblate shapes. For comparison, 
the densities at the spherical shape are also displayed. It 
can be seen that significant central depressions or central 
semi-bubble appear in the spherical case. However, the 
situation is considerably altered even with a small shape 
change. Only weak central depressions are seen at the 
small deformations, see Fig. 10. To have a further under- 
standing, we calculated the corresponding single-particle 
potentials in the z-axis, shown inside Fig. 10. It can be 
seen that the proton potential has a considerable change 
with changing the deformation. The spherical proton po- 
tential has a significant hump at the center of the nucleus. 
This implies that the Coulomb energy can be consider- 
ably reduced by forming the center semi-bubble at the 
spherical shape. The Coulomb energy can also be re- 
duced by generating the deformation of the nucleus. In 
the reduction of the total energy of the nucleus, there 
is competition between forming the central semi-bubble 
and generating the deformation. The deformation can 
affect the shell structure and then density distributions, 
and vice versa. In a self-consistent model, such as the 
SHF approach, densities are fed back into the potential, 
which amplifies the coupling between deformations and 
densities. 



III. SUMMARY 

In summary, the density distributions of superheavy 
nuclei have been investigated with the Skyrmc-Hartree- 
Fock model. To test the model and parameters, the a- 
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FIG. 6: Calculated density distributions of nuclei of the deformed 270 > 280 ' 290 no. The left and middle columns show the 
distribution profiles of p(z, 0) along the z-axis and p(0, r) along the r-axis, respectively. For comparison, distributions at 
assumed spherical shape are also displayed in the right column. 



decay energies of even-even superheavy nuclei and the 
charge density of 208 Pb are calculated and compared to 
existing experimental data. For axially-symmetrically 
deformed nuclei, the density distributions in the sym- 
metric z- and the vertical r-axes or the two-dimension 
distributions were calculated. The distribution in differ- 




FIG. 7: Two-dimensional proton density contours of the nu- 
cleus 280 110, calculated with the SkI4 force. The numbers on 
contours are the values of densities in fm -3 . 



ent directions can be different. The high-j low-j z orbits 
have important contributions to the densities at nuclear 
surface in the z-axis, while high-j high-j z orbits have 
important contributions at surfaces in the r-axis. The 
deformation effect was found to be significant in the cal- 
culation of the density distribution in the 292 120 nucleus. 
Only a weak central depression was seen in the deformed 
case for 292 120, compared to the predicted semi-bubble 
at the spherical shape. 
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FIG. 8: The SkI4-calculated contributions (in z-axis) of the 
j*=l/2 orbits in the 280 110 case, with the modifications of 
pairing occupation probabilities. The dashed line represents 
the summed contributions of the orbits. 
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FIG. 9: Calculated energy curves of 292 120 with different 
Skyrme forces, as a function of quadrupole deformation. 
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FIG. 10: The SkI4 calculated density distributions p(z,Q) and p(0,r) in 292 120 at f3 2 = -0.12 (a), 0.11 (b) and 0.0 (c), 
corresponding to the energy minima shown in Fig. 9 and the spherical shape. The inset figures show the single-particle potentials 
in the z-axis. 



